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Abstract 


A new analytic result in acoustics called “Formulation IB,” proposed by Farassat, is used to 
compute the loading noise from an unsteady surface pressure distribution on a thin airfoil in the 
time domain. This formulation is a new solution of the Ffowcs Williams-Hawkings equation with 
the loading source term. The formulation contains a far field surface integral that depends on 
the time derivative and the surface gradient of the pressure on the airfoil, as well as a contour 
integral on the boundary of the airfoil surface. As a first test case, the new formulation is used 
to compute the noise radiated from a flat plate, moving through a sinusoidal gust of constant 
frequency. The unsteady surface pressure for this test case is analytically specified from a result 
based on linear airfoil theory. This test case is used to examine the velocity scaling properties 
of Formulation IB and to demonstrate its equivalence to Formulation 1A of Farassat. The new 
acoustic formulation, again with an analytic surface pressure, is then used to predict broadband 
noise radiated from an airfoil immersed in homogeneous, isotropic turbulence. The results are 
compared with experimental data previously reported by Paterson and Amiet. Good agreement 
between predictions and measurements is obtained. Finally, an alternative form of Formulation IB 
is described for statistical analysis of broadband noise. 
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Nomenclature 


A c = correlation area (m 2 ) 
b = airfoil semi-span (m) 

C = airfoil chord (m) 

Co = ambient sound speed (m/s) 

/ = frequency (Hz) 

/ = geometry function for airfoil surface (Fig. 1) 

E* = complex combination of Fresnel integrals (Eq. (9g)) 

g = velocity-to-pressure transfer function 

1C = time independent kernel function (Eq. (24) ) 

k = u/U, convective wave number (m _1 ) 

k = kCj 2, reduced frequency 

C\ — streamwise integral length scale (m) 

L c = correlation length (m) 

M = F/co, Mach number vector 
M r — M • r/r, Mach number in radiation direction 
M v — M • Mach number in direction of v 
A P— unsteady airfoil surface pressure jump (Pa) 
p = unsteady airfoil surface pressure (Pa) 
p' = sound pressure radiated to observer (Pa) 
dp/ds = surface pressure gradient in the direction of surface velocity (Pa/m) 

R = length variable associated with radiation distance in retarded time (Eqs. (23b) and (23c)) 
R ww = correlation function of w (Eq. (19g)) 
r = x — y , sound radiation vector (m) 

Sww — power spectral density of w in ft-space (m 3 /rad-s 2 ) 
t = observer time (s) 

U = uniform freest ream speed (m/s) 
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u = unsteady streamwise velocity (m/s) 

V = airfoil velocity vector in reference frame fixed to undisturbed medium (m/s) 
w = unsteady upwash velocity (m/s) 

x = [x\,X 2 , # 3 ] T , observer position in reference frame fixed to undisturbed medium (Fig. 1), 
fixed to airfoil reference frame in Section 5 only 

V — [yi> 2 / 2 , 0] T , surface source position in reference frame fixed to undisturbed medium (Fig. 1), 

fixed to airfoil reference frame in Section 5 only 

p = Vl-M 2 

A = Co//, acoustic wavelength (m) 
p = Mk/p 2 

'ip = directivity angle (Fig. 7) 

v — unit inward facing normal on surface edge (Fig. 1) 

9 = angle between airfoil surface normal and radiation direction (Fig. 1) 
po = ambient density (kg/m 3 ) 
r = source time (s) 

$ = random phase variable 

U3 = 27 r/, circular frequency (rad/s) 

Subscripts 

1,2,3 = chord wise, spanwise, and normal coordinate directions (Fig. 3) 

1A = solution calculated by Formulation 1A 
IB = solution calculated by Formulation IB 
ret = evaluated at retarded time t — v/cq 


4 of 51 



1- Introduction 


The development of analytical methods to predict noise radiated from an airflow over a rigid 
body has been an active research topic within the aeroacoustic community for decades. Research in 
this area has, in large part, been motivated by the desire to incorporate the results of aeroacoustic 
analysis into an aerodynamic design methodology. The present work is similarly motivated, and 
the resulting formulation should lend itself well to an engineering design tool suite. 

The current work is specifically focused on the calculation of far field noise that results from 
fluctuating pressure on a solid surface. The acoustic analogy [1] provides a framework for developing 
methods to predict noise from many types of sources, including noise due to unsteady surface 
loading. Such noise is mathematically described by the loading source term, or “dipole term,” of 
the Ffowcs Williams-Hawkings (FW-H) equation [2]. Because the noise due to an airflow over a 
rigid surface is typically dominated by dipole radiation, the acoustic formulations of interest in 
this work are determined by solutions of the FW-H equation with the loading source term, i.e. 
neglecting the thickness and volumetric terms. 

The solution of the FW-H equation can be written in many forms. A new solution, called 
“Formulation IB,” is presented herein. Most types of broadband noise, including noise that is 
dominated by an airfoil’s leading and trailing edges, can be calculated with the proposed Formu- 
lation IB, which is the simplest loading noise prediction formula known to date. This simplicity 
makes the new formulation highly suitable for statistical analysis of broadband noise for rotating 
surfaces. 

In Section 2, Formulation IB is derived for the specific case of a flat plate undergoing a general 
nonuniform motion. For low Mach numbers and distant observers, the dominant term in this 
formulation is a far field surface integral that depends on the time derivative and the gradient of 
the airfoil surface pressure. The formulation also contains a contour integral on the boundary of the 
airfoil surface. This line integral vanishes along the trailing edge if the Kutta condition is imposed. 

In Section 3, Formulation IB is used to calculate the noise radiated from a flat plate moving 
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through a sinusoidal gust of constant frequency. The unsteady surface pressure used in this test 
case is an analytical result from linearized airfoil theory that is taken from the work of Amiet 
[3, 4]. A mesh refinement study is performed to demonstrate the equivalence of Formulation IB 
with Formulation 1A [5], a previously developed acoustic formulation that is also a solution of 
the FW-H equation. Results from this analytical test case are also used to examine the velocity 
scaling properties of Formulation IB, which are found to be consistent with the results of Curie [6] 
and Ffowcs Williams and Hall [7]. The directivity of the noise induced by a periodic gust is also 
examined. 

In Section 4, the single- frequency surface pressure in Section 3 is extended by spectral rep- 
resentation to serve as an analytic broadband source model for incident turbulence noise. This 
stochastically modeled surface pressure is used in Formulation IB to predict broadband noise 
caused by the interaction of a NACA 0012 airfoil and homogeneous turbulence. The resulting 
calculations are compared with experimental data previously reported by Paterson and Amiet [8]. 
In Section 5, an alternative acoustic formulation is described for statistical analysis of broadband 
noise. 


2. Acoustic Formulation 

Consider a flat, finite surface moving in the plane x$ = 0 along a velocity vector V. The velocity 
vector V and the surface geometry are related to the coordinate axes as pictured in Fig. 1. Let 
f(x i,#2, t) denote a geometric function where / = 0 on the surface edge and / > 0 on the interior 
of the surface. Let v = V / denote the unit normal which lies in the plane of the surface, is normal 
to the edge, and is directed inward (Fig. 1). Note that V need not be constant in space or time. 
The only stipulation on the velocity is that the motion of the surface is in the same plane as the 
surface. 

Let x = [x\^X 2 - ) x^\ r denote the position of an observer, and y = [yi,2/2? 0] T denote the position 
of a source point on the plate’s surface (Fig. 1). The spatial frame of reference is fixed to the 
undisturbed medium. The unsteady surface pressure p(y, r) gives rise to sound that initiates at 
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source time r and radiates to the observer along f = x — y. This sound is described by p f (x, £), 
the perturbation pressure that arrives at the point {xi,X 2 : xs) at time t. For some of the subtle 
mathematical details that are involved in the following derivation, see [ 9 ]. 

From the FW-H equation, the loading noise is given by a solution of 

^ ^ - v 2 p' = - v • [p n if (/) * (ar 3 ) ] (la) 

where Co is the ambient sound speed and n is the unit surface normal that, for the present case, is 
equivalent to £3, the unit vector in the direction of the ^3-axis. The function H(f) is the Heaviside 
step function and 8(x 3) is the Dirac delta function. Evaluating the divergence in Eq. (la) yields 

4 - v V = -p(xi,x 2 ,t)H(f)S'(xs) (lb) 

Cq ot z 

where 8 l (x$) denotes differentiation with respect to #3. 

Eq. (lb) is the wave equation with an inhomogeneous source term, and its formal solution in 
an unbounded domain requires the Green’s function <5(#)/47rr, where g — t — t + r/co , and t and 
r are the observer and source times, respectively. The solution of Eq. (lb) can then be written in 
the form 

47T p'(x,t) = -f [ ^-p{yuy2,r)H{f)5'{yz)dydT 

J —00 ^ 

Let t — y g and integrate with respect to g . The result can be written 

47T p{x,t) = -[ - [p] ret H(F) 5'{yz) dy ( 2 a) 

where the subscript “ret” denotes evaluation at retarded time t — r/co, and F is 

F(yi,y2 ; x,t) = f(yi,y2,t- r/c 0 ) = [/] re t (2b) 

Integration with respect to on the right-hand side of Eq. ( 2 a) yields 

47 T p'{x,t) = f 4 “ { H (F) | %1 d V2 ( 3 a) 

JW %3 l r ) ys=0 
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Performing the differentiation in the integrand of Eq. (3a) yields 


^ J IZk (jT-) j = —[f 3 p] ret H(F) + \[hpUH(F) + ~ 

oyz \ r J cor r 1 cor 


df - 

fr r3P 


5(F) (3b) 


where p is the time derivative of pressure evaluated by an observer that is fixed with respect to 
the medium at rest, and r% is the third component of the unit radiation vector r = r/r. Clearly, 
then, ^3 = es * f = cos 6 , where 6 is the angle subtended by the surface normal and the radiation 
vector, as shown in Fig* 1. The first and second terms on the right-hand side of Eq. (3b) are of the 
similar form QH(F). In [10], it is shown that the integration of these two terms in Eq. (3a) can 
be written 


[ QH(F)d yi dy 2 = f QdZ= f % dS (4a) 

JSft 2 JF > 0 Jf > 0 l 1 — M r Jret 

where eE is the surface area element of the acoustic planform of F > 0 . Also, M r — M * r is the 
Mach number in the radiation direction, where M = V/cq is the local Mach number vector of the 
surface. 

The integrated value of the third term on the right-hand side of Eq. (3b) is determined as 
follows. This integral is of the form 


1 = 



q(yuy2)5(F)dyidy 2 


(4b) 


The differential surface element dy\dy 2 can be written [10] 


dyi dy 2 = dC dAf = 


dCdF 

|V 2 F| 


(4c) 


where dC and dAf are differential elements of arclength that are, respectively, parallel and normal 
to the the surface edge defined by F = 0 as shown in Fig. 2. The notation V 2 denotes the surface 
gradient in the y\ 2 / 2 -phtne. Moreover, it can be shown [10] that 


dC _ d£ 

V 2 -F | [ 1 — Af r ] re t 


(4d) 


where d£ is an element of arclength along the surface edge defined by / = 0 . Eq. (4b) can then be 
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written 


= [ q(yu 2 / 2 ) S(F) 

JSR 2 


dCdF 


diVuVi) 


dC = 


qiyum) 


di 


(4e) 


| V 2 F | Jf= 0 I V 2 F I .//= 0 [ 1 — Jret 

Note that the surface time derivative df /dr in Eq. (3b), and contained in < 7 ( 2 / 1 , 2 / 2 ) in Eq. (4e), 
is referenced to the undisturbed medium. However, df /dr can be related to the local normal 
in-plane velocity of the edge contour as 


D£ 

Dt 


f = 0 


= ^- + V -Wf =^- + V -0 = 0 
or or 


(5a) 


where the subscript “/ = 0” denotes evaluation on the edge contour. This relation yields df /dr = 

-9 -P. 


All three terms on the right-hand side of Eq. (3b) are now integrated in Eq. (3a) over the 
physical surface f(xi^X 2 : t) > 0 by using Eqs. (4a) and (4e). Note that in Eq. (la), p is referenced 
to the medium at rest, e.g., as measured by a transducer that remains stationary as the surface 
passes by it. The quantity p can also be related to dp /dr , the time derivative of pressure in the 
reference frame of the moving surface, e.g., as measured by a transducer attached to the surface. 
This relation is 


• dp dp 

p= ^~ v ¥. 


(5b) 


where dp/ds is the gradient of p in the direction of V, V is the local magnitude of V, and s is in the 
direction of the velocity V of the surface in the reference frame fixed to the undisturbed medium. 


Incorporating all of the above results into Eqs. (3a) and (3b), the solution of Eq. (la) can be 
written 


Attp(x : t) = 


f > o L 


( dp/ dr — V dp/ds ) cos 0 


+ 


/>o L 


c 0 r ( 1 - M r ) 
p cos 9 


dS 

J ret 


(6) 


r 2 ( 1 - M. 


w J ret 


dS- 


f = o L 


M v p cos 0 
r (1 — M, 


w J ret 


di. 


where M v = the Mach number in the direction of is. Eq. ( 6 ) will be referred to as Formulation 
IB. 
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The first and third integrals in Eq. (6) are far field terms, and the second integral is the near field 
term. The relative contributions of the terms in Eq. (6) are now considered, under the conditions 
of low Mach number and an observer in the acoustic far field, i.e., 

M < 1 , r > A (7) 

where A is a typical acoustic wavelength of interest. With respect to M and r, the surface far 
field integral, i.e., the first integral in Eq. (6), is proportional to 1/r, whereas the second and third 
integrals are proportional to 1/r 2 and M/r, respectively. Therefore, the far field surface integral 
dominates the signal under the conditions in Eq. (7). 

Again, Eq. (6) is valid for the case of nonuniform surface motion. Therefore, Eq. (6) can be 
used, as is, to predict loading noise from rotating surfaces. Its predecessor, Formulation 1A [5], 
is significantly more complicated in its rotational form and cannot be approximated by only one 
surface integral in the far field. Such a significant simplification for far field calculations makes 
Formulation IB more suitable for statistical analysis of broadband sources for rotating surfaces. A 
statistical formulation based on Eq. (6) will be addressed in Section 5. However, the focus of the 
current work is the application of Formulation IB in the form of Eq. (6), as will be demonstrated 
in the following two sections. 

3. Sinusoidal Gust of Constant Fre q uenc y 

Any noise prediction made with Eq. (6) will be only as good as the input surface pressure 
p(y, r). Such time-dependent pressure data could be taken from experimental measurements or 
a computational fluid dynamics (CFD) calculation. However, in this initial work, an analytic 
expression is used for p(y, r) to serve as input data for the new acoustic formulation. To this end, 
an analytic formulation from thin airfoil theory is used to describe the unsteady surface pressure that 
results from a sinusoidal gust of constant frequency. This particular surface pressure formulation is 
chosen as an example to establish the equivalence of Formulation IB and Formulation 1 A [5] and to 
examine the velocity scaling properties of Formulation IB. In addition, the frequency dependence 
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of the far field directivity is discussed. 


3.1 Surface Pressure from Thin Airfoil Theor y 

Consider a rectangular flat plate, in rectilinear motion, as in Pig. 3. The velocity vector V = 
[— C/,0, 0] T , where U is a constant subsonic speed. For the following examples, the plate’s surface 
and its boundary, f > 0, are defined by the rectangle { 0 < aq < C } x {— b < X 2 < b} in the 
plane xs = 0 . This surface is assumed to have an unsteady pressure distribution that is analytically 
prescribed from thin airfoil theory, as discussed below. 

Amiet [3, 4] presents closed form expressions for the unsteady pressure on the surface of a 
thin airfoil of infinite span. The airfoil is assumed to move rectilinearly through a sinusoidal 
gust. Analytical methods are used to solve the two-dimensional, time-dependent linear potential 
equation by representing the solution as a product of spatial and temporal solutions. The solution 
is represented as a truncated series in which higher order terms are neglected. 

A complex-valued representation for the airfoil surface pressure is assumed to arise from a 
stationary gust in one spatial dimension. This gust can be written in the stationary variable 
x\ — Ut as 

w(x 1 - Ut) = wo (8) 

where k = uj/U is the streamwise convective wave number, and wq is the gust amplitude. This 
gust and the airfoil surface pressure that it generates are, for now, considered functions of a single 
amplitude and frequency. 

The unsteady pressure that arises at a point (yi, 2/2,0) on the airfoil surface as a result of the 
gust in Eq. (8) can be written 


AP(yi,f) = 2p 0 Uwog(yi,k)e tkm (9a) 

where po is the ambient density, k = kC/2 is the reduced frequency (based on the semi-chord), and 
g(y i, k) is a transfer function whose form is dependent on the frequency of interest. The factor of 
two in Eq. (9a) indicates that the pressure is assumed to be antisymmetric between the upper and 
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lower surfaces, and this expression thereby accounts for the pressure on both sides of the airfoil, 
i.e., the pressure jump. In [8], the suggested parameter to delineate between the low and high 
frequency regimes is (i — Mk/f3 2 , where /3 = \/\ — M 2 , with a demarcation value of fj, = 0.4. 

For low frequencies, fi < 0.4, the transfer function is 


g(y uk) 


l 


yi 


G s {k*y h ' q{yuM) , 


H < 0.4 


(9b) 


where k* = 
mated by 


A;//? 2 , Gs is the classical Sears function [11], which, for the present work, is approxi- 


G s m 


i 

1 + 2.4&* 


1 

2 


+ 2 irk* 


(9c) 


as suggested in [8], and 


q(y u M ) = M 2 (2j/i/C- 1) + (1-0) In M + /3 ln(l+0) - In 2 (9d) 


For high frequencies, /i > 0.4, the transfer function is the sum of a leading edge solution and a 
trailing edge correction [4], i.e., 


g(yi ,k) = (gi+ 9 2 ) e~* , „ > 0.4 


(9e) 


where 


gi (yi,k) 


1 

[27rfcj/i(l + M)/C]z 


92(yi,k) 


— 1 + ( 1+0 E*[4:fi(l — 3/1/C) ] 
[27Tfc(l + M)]5 


(9f) 


and 

E*(0 = f — — — du = C(0 - i S(0 (9g) 

Jo ( 2 nu )2 

The quantities C(£) and <S(£) are the Fresnel cosine and sine integrals and will be numerically 
evaluated by the formulas in [12]. The final representation for the unsteady surface pressure p(yi, r), 
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assumed to be a real quantity, is 


p(yi,r) = U{-AP(yi,r)} (9h) 

The pressure jump is negative in Eq. (9h) because the acoustic formulation was derived from 
Eq. (la), in which the unit surface normal £3 is assumed to point into the fluid, i.e., in the positive 
X 3 direction on the upper surface, and in the negative x% direction on the lower surface. Therefore, 
using the same positive surface normal on both sides of the airfoil, the sum of the pressure on both 
sides is A P = _P upper — P| ower , and this expression is the negative of the conventional notion of a 
pressure jump. 

Some differences are noted here between Eqs. (9e) and (9f) and their counterparts in [3] and 
[8]. Such differences include the choice of the coordinate-axes origin and the spatial normaliza- 
tion employed by the authors. Another difference arises from the use of complex conjugates in 
the present paper, which serves to make the present notation more consistent with the eventual 
broadband formulation in [8], which is employed by the present authors in Section 4. The above 
surface pressure was originally proposed in order to derive an expression for unsteady lift that was 
ultimately incorporated into a frequency domain acoustic formulation [3, 8, 13]. However, in the 
present work, the unsteady pressure jump itself will be used as input to Eq. (6) for a time domain 
prediction. Finally, note that for both transfer functions, the pressure jump vanishes at the trailing 
edge x — C\ i.e., the Kutta condition is satisfied. The correctness of this trailing edge boundary 
condition has been a subject of debate for many years, for example as reviewed by Howe in [14]. 
However, as the response of the airfoil leading edge is expected to dominate the noise that arises 
from an incident gust, determining the most appropriate treatment at the trailing edge is not an 
issue for the current problem. 

3.2 Grid Refinement Stud y 

The surface pressure in Eqs. (9a)-(9h) is now used to numerically demonstrate the equivalence 
of Formulation IB and and the loading noise terms of Formulation 1A [5]. Formulation 1A forms 
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the basis of WOPWOP, a rotor noise prediction code developed at NASA Langley Research Center 
[15]. To predict loading noise from an airfoil in uniform rectilinear motion, Formulation 1A is 
simplified to 



dp /dr cos 0 
c 0 r(l-M r )\ 


dS 

ret 



P COS0 

_r 2 ( 1 — M r ) 2 _ 


dS- 

ret 



p cos0 (M r — M 2 )! 

| r 2 ( 1 — M r ) 3 J ret 


( 10 ) 


The entire formulation here is integrated on the surface interior. At first glance, the form of Eq. 
(10) appears no more complex than Eq. (6), but only because of the simplicity of uniform, rectilinear 
motion. Unlike Eq. (6), Eq. (10) is not applicable to a rotating surface. The loading noise terms of 
the full Formulation 1A are indeed applicable to rotational flow, but the full formulation is more 
complex than Eq. (10) (see [15]). 

The far field noise radiated from a thin airfoil in a one- dimensional, single-frequency sinusoidal 
gust is now calculated with Eq.s (6) and (10). Let p\ B and p f 1A denote the sound calculated by 
Eqs. (6) and (10), respectively. If the input surface pressure p(y, r) is exactly specified at every 
point on the airfoil surface, then the only error in the numerical solution of Eqs. (6) and (10) arises 
from the quadrature formula that is chosen to perform the surface and contour integrations. In this 
case, the midpoint rule is the quadrature of choice. Therefore, the equivalence of Eqs. (6) and (10) 
is demonstrated if the difference \p\ B — p\a I diminishes in mesh refinement like the cumulative 
error expected from the midpoint rule, i.e. that the error is proportional to the square of the mesh 


spacing. 


The plate’s rectangular dimensions are determined by a chord length of C = 0.5 m and a span 
of 26 = 2.0 m. The plate is moving at a Mach number of 0.2, and the sound speed is taken to 
be 343 m/s. The ambient density is po = 1.23 kg/m 3 , and the upwash amplitude is wo = 0.05 U. 
The observer position for this test case is x = [—1, 0, 1] T , in meters. Fig. 3 roughly depicts the 
direction of this observer position, although the distance is not to scale. 


The calculation is performed for one time period of the surface pressure fluctuation at frequen- 
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cies of / = 25 Hz and / = 1 kHz, with 32 time steps in each period. These choices of frequency, at 
the prescribed observer location, will test both the near field and far field equivalence of the two 
formulations. Note that the transfer functions in Eqs. (9b) and (9f) are singular at y\ = 0, and 
the spatial derivative of Eq. (9f) is singular at y\ — C . Although both singularities are integrable 
in an analytic sense, they would cause the quadrature error to deviate from that of the midpoint 
rule which, by its definition, requires sufficient smoothness throughout the interval of integration. 
Therefore, the domain of streamwise integration is restricted to an interval of the form 

eC < y\ < (1 — e) C (11) 

where e is a small, positive parameter. 

Each calculation is performed on a sequence of six surface grids: {10 x 40}, {20 x 80}, {40 x 
160}, {80 x 320}, {160 x 640}, and {320 x 1280}. Because the surface pressure is cast in only one 
spatial variable yi, and the observer location is symmetric relative to the airfoil span, the acoustic 
prediction is relatively insensitive to the discretization in ?/2- These discretizations are employed 
as a more general example of the types of resolution that would be required for a case in which 
spanwise pressure variation would significantly affect the far field acoustics. 

The grid is clustered via Roberts’ transformation [16] near the leading and trailing edges of the 
plate in order to accommodate the parameter e = 0.02 on the coarser meshes. The 20 x 80 grid 
is shown in Fig. 4. The maximum values of \p\ B — p\^ | during each time period are shown as a 
function of the number of grid points on a log-log plot in Fig. 5. The abscissa Nc is the number of 
surface elements in the chordwise direction. The dashed line represents a fictitious quantity whose 
values are specifically calculated to be directly proportional to TV^ 2 . The slopes of both calculations 
are visibly parallel to a slope of -2, thereby demonstrating that 


p'ib ~P'iA I = 0(N C 2 ) = 0(Ay?) 


which is the cumulative error expected from the midpoint rule of integration. 
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3.3 Velocity Scalin g Laws 


Attention is now turned to the way in which the intensity of the far field noise, as determined by 
Eq. (6), will scale as a function of velocity when the surface pressure is described by Eqs. (9a)-(9h). 
The conditions of a far field observer in a low Mach number flow (Eq. (7)) will be assumed. The 
Mach number range of interest is 0.01 < M < 0.2. Velocity scaling laws will be determined for 
both compact and noncompact sources. 

First, the source is assumed to be compact, i.e., C < A. This condition will be achieved by 
requiring the frequency to be proportional to velocity, / ~ U/C , for a sufficiently low range of 
frequency. This assumption is consistent with the compact source argument put forth by Curie [6] . 
If a frequency range of 4 Hz < / < 80 Hz is chosen for the Mach number range 0.01 < M < 0.2, 
then the acoustic wavelength A ranges from 85.75 m at M = 0.01 to 4.29 m at M = 0.2, for 
an ambient sound speed of 343 m/s. If the plate’s physical dimensions are as in the above mesh 
refinement problem ( C = 0.5 m), then the ratio \/C ranges from 8.58 at / = 80 Hz to 171.5 at / = 
4 Hz. These ratios are assumed sufficient for an acoustic source to be considered compact. Note 
that (i < 0.4 throughout these ranges of Mach number and frequency. 

The calculation is performed on a 100 x 400 surface grid, with the streamwise integration interval 
restricted as in Eq. (11), with sufficient grid clustering near the leading and trailing edges to allow 
for € = 0.003. The observer position is chosen at a distance of 100 m, directly above the plate’s 
center, i.e., x = [0.25,0, 100] T in meters. This location places the observer in the acoustic far field 
for the entire range of frequency. 

The upwash amplitude is wo = 0.05 U for each of 50 equally spaced Mach numbers between 0.01 
and 0.2. A separate calculation is run for each Mach number and its corresponding frequency. The 
surface pressure in Eqs. (9a)-(9d) and (9h) is used as input to Eq. (6) to predict the far field sound 
p f (x,t). Each calculation is performed for one acoustic period T of the corresponding frequency, 
with 64 time-steps. The average intensity I(x) of the acoustic signal at the observer x , assuming 
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spherical spreading, is then calculated by 

/(f) = - f T [p'(^)] 2 ^ 

T Jo Po co 

The average acoustic intensities for a compact source, as a function of Mach number, are repre- 
sented as squares in Fig. 6. The slope of these results on a log-log plot can be visually determined 
by observing their proximity to the dashed line with the slope of six. This U 6 proportionality is 
consistent with Curie’s result [6], as expected from the conditions placed on the calculations. 

The demonstration of a velocity scaling law is next desired for a noncompact source where 
H > 0.4. Therefore, the restriction that / ~ U/C must be lifted so that / is independent of U. 
The simplest such condition is that / is constant, in which case the Strouhal number remains a 
function of U . A series of calculations is again performed, as above, with the only parameter change 
being that the frequency is held constant at 1 kHz throughout the range of Mach number. The 
input surface pressure for this case is given by Eqs. (9a) and (9e)-(9h). The computed acoustic 
intensities at 100 m are represented by circles in Fig. 6. In this case, the acoustic intensity scales 
approximately as £/ 5 , a result that is consistent with Ffowcs Williams and Hall [7]. 

3.4 Directivit y 

As a final exercise in this section, the directivity of a single- frequency source is examined. The 
radiated noise p f (x : t) is calculated at many locations on a circular arc in the plane X 2 = 0 that is 
geometrically centered on the plate’s upper surface, as shown in Fig. 7. The arc trajectory (r, ip) 
is determined by r — 3 m and 0 < ip < tt. The directivity is determined by the peak pressure 
amplitude |p'| max that is calculated at each position on the circular arc during one period in time 
for a given frequency. The flat plate’s dimensions and surface discretization are the same as in 
the previous example, and there are 128 time-steps in a period. The observer path, 0 < ip < tt, 
is discretized into 128 equally spaced locations. The freestream Mach number is 0.2 and the gust 
amplitude is wq = 0.05 U. Fig. 8 shows the results, in polar form, for frequencies of 1 kHz and 2 
kHz. The higher frequency results are scaled to show both loci on the same plot. As expected, the 
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2-kHz solution has twice as many lobes as the 1-kHz solution. Furthermore, the general downstream 
directivity of the lobes is consistent with the results of previous authors e.g., [17, 18]. Finally, note 
the frequency dependent positions of the lobes with respect to a fixed observer. For example, an 
observer at ^ = 90° will receive a signal that is near the peak amplitude for the central lobe of 
the 1-kHz signal, whereas the 2-kHz signal is near a local minimum for the same observer. This 
frequency dependent characteristic of directivity is discussed further in the next section. 

4. Broadband Prediction with Com p arison to Ex p eriment 

The analytic surface pressure in the previous section is extended to model a broadband source 
on a slender airfoil at zero angle of attack. This broadband surface pressure is used as input 
to Formulation IB to predict far field radiation, and the results are compared with experimental 
measurements. Amiet [19] has previously proposed a broadband solution to this problem in the 
time domain by using the transfer functions in Eqs. (9e)-(9g). His resulting solution was a Fourier 
transform of the frequency domain solution in [8] and [13]. In the current work, the high frequency 
formulation in Section 3 is explicitly extended to a broadband source application with user- specified 
spectral content. 

4.1 Ex p eriment Descri p tion 

The experiment that is modeled in this section is reported by Paterson and Amiet [8]. A NACA 
0012 airfoil is placed between two vertical plates at zero angle of attack in the test section of an 
open jet wind tunnel. A schematic of this experimental setup is shown in Fig. 9. The airfoil has a 
chord length of 0.23 m and a span of 0.53 m. Turbulence is generated by a grid upstream of the 
airfoil. Noise propagates from the test section into an anechoic chamber that is instrumented with 
six microphones. The microphones are located on the tunnel centerline, on an arc of radius 2.25 m, 
relative to the airfoil’s geometric center. The microphone location of interest here is at 90 degrees 
(Fig. 9). Far field noise measurements of the incident turbulence on the airfoil are determined by 
subtracting microphone measurements, with and without the model in the test section, at each of 
five tunnel speeds: 40, 60, 90, 120, and 165 m/s. 
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4.2 Broadband Anal y sis 

For prediction purposes, the airfoil is modeled as a flat plate in a periodic gust which gives 
rise to an unsteady surface pressure that is a broadband extension of the analytic formulation in 
Section 3. The airfoil geometry is oriented with respect to the coordinate axes as in Fig. 3, with 
{0 < xi < C} x { — b < X2 < b}, where C = 0.23 m and 2b = 0.53 m. The components fci, k 2, 
and &3 of the wave number vector are oriented along the axes xi,X2, and x$ , respectively. As 
encountered by the airfoil surface in the plane x$ = 0, the normal component of the turbulent 
velocity field can be written 

/ oo r oo 

/ w{k 1 M)e- i[kl{xi - m)+k2X2] dk 1 dk 2 ( 12 ) 

-OO’/ — OO 

where w(k \ , fe) is the wave number spectrum of the gust amplitude, defined as the inverse Fourier 
transform of w{x\, x^, t). Furthermore, this two-component velocity spectrum is interpreted as a 
fully three-dimensional spectrum whose form results from the integration of all wave numbers that 
are normal to the airfoil surface: 

/ oo 

w(kiM,h)dh 

-OO 

The comp lex- valued, unsteady surface pressure jump that arises at a point (yi, 2/2,0) on the airfoil 
surface from the incidence of a turbulent velocity field of the form in Eq. (12) is given by 

/ oo roc 

/ w{kiM)g{y\MM)e l(<klUt ~ k2y2) dkidk 2 (13) 

-oo J — oo 

A straightforward approach for predicting the desired broadband far field measurements is to 
use the real part of Eq. (13) as input to Formulation IB. The velocity spectrum required in this 
approach can be derived from well known power spectrum formulas for homogeneous, isotropic 
turbulence, such as von Karman’s formula [20]. A formulation for the dual wave number transfer 
function giyukifo) is given by Adamczyk [21]. 

However, because one of the objectives of the current work is to reproduce the results of Paterson 
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and Amiet [8] from a time domain perspective, an approach similar to that taken in [13] will be 
taken to model the surface pressure. Amiet [13] argues that, within certain limitations, integration 
over all spanwise wave numbers is not required. His conclusion, derived mathematically in the 
frequency domain, is that only one spanwise wave number contributes to the sound detected by an 
observer in a given location. In particular, Amiet focuses on an observer in a spanwise symmetric 
location, for which only the zero spanwise wave number needs to be considered. This result is 
argued to be exact in the limit of infinite span and a good approximation for an airfoil of finite 
span that encounters a high frequency disturbance. 

The derivation of Amiet’s analytical result can be generally described as follows. First, Eq. 
(13) is transformed into Fourier space, with the stationary variable x\ — Ut being replaced by 
oj = kiU. Then, a two-point cross-correlation function is formed and related to the far field power 
spectrum through Kirchhoff’s formula [22] and Curie’s result [6]. In order to follow a similar line 
of reasoning in the time domain, Eq. (13) itself must be related to the far field acoustic pressure 
through Formulation IB. In the case of a distant observer directly overhead of a finite-span airfoil, 
the terms r, M r , and 0 in Eq. (6) are weak functions of yi and 7/2? and therefore will be considered 
constants. Furthermore, for the observer position considered here, the differences in retarded time, 
as a function of airfoil surface location, can be neglected. These assumptions are consistent with the 
acoustic model employed by Amiet [13]. For the present problem, including the above assumptions, 
Eq. (6) is approximated by 


4717 /( 0 ?, t) 


cos 9 


cor (1 - M r ) 

M cos# rb 


r(l-Mr) J- b 


rrd _ d 

Jo J b dt P ^ VuV2 ' 1 ~ ry/C °^ + U dyl P ^ Vuy2 ' 1 ~ dy 2 dyi 

p(0,y2,t-r/co)dy 2 (14) 


where the over-bars on 0, r, and M r denote mean values over the airfoil surface, and therefore 
the retarded time t — f/co is constant for fixed t. The near field term (second integral in Eq. (6)) 
has been deleted because of its 1/r 2 proportionality. However, the third integral must be retained 
because the Mach number is not necessarily assumed to be small. Only the leading edge line integral 
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survives the contour integration in the third term in Eq. (14) because M v — 0 along the airfoil tips 
and the Kutta condition is assumed along the trailing edge. 

Before the surface pressure p is specified, Eq. (14) is further simplified. For convenience, the 
terms 0 and M r will be neglected, as they are small (M r « 0 and cos0 « 1) for a distant observer 
directly above the airfoil. With these additional simplifications, if —A P in Eq. (13) is substituted 
for p(y\, t/2, t ~ r/co) in Eq. (14), the far field acoustic pressure can be approximated in the form 


4717 / (x, t) 


2p 0 U 


C pb poo poo 


+ 


/oo 

T\ ( 2/1 ,k\,k 2 ) e^ fcl ut ~ k ^) dk 2 dki dy 2 dy\ (15a) 

-OO 

2p0 UM f h 7°° 1 JF 2 (0, k \ , k 2 ) dk 2 dk\ dy 2 
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where 


Fi(yi,h,k 2 ) 

!F 2 (0, ki,k 2 ) 


-w(ki,k 2 ) 


ih U g(yi,ki,k 2 ) + U 


iMM) 


e ~ikir/co 


-w(ki,k 2 )g(0,h,k 2 )e * feir/c ° 


(15b) 


Sufficient conditions [23] on A P and its derivatives have been assumed for the commutation of 
integration and differentiation in Eqs. (15a) and (15b). The y 2 integration in Eq. (15a) can now be 
explicitly evaluated, yielding 


47T p f (x : t ) 
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C poo poo 
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cor 
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2sin(Ar 2 5) 

k 2 

2 sin(k 2 b) 

k2 


Fi(yi, k\,k 2 ) e lklUt dk 2 dk\ dyi 


JF 2 (0, fci,fc 2 )e* felC/t djt 2 dfci 


(16) 


Integrating with respect to fe, the term sin(&2&)/fc2 acts like a Dirac delta function when integrating 
over an unbounded domain, and the result is 


47T p'(x,t) 
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/ 2tt F\( yi, ki,0) e iklUt dk\ dy\ 

cor J 0 J 

r 

— OO 

1 2poUM 
f 

POO 

/ 2n T 2 (0, k\ , 0) e i klUt dk\ 

J — 00 


(17) 
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Eq. (17) indicates that only the zero spanwise wave number contributes to the noise detected by 
the far field observer. Eqs. (16) and (17) are time domain analogies to Eqs. (15) and (17) in [13]. 
Furthermore, Eq. (17) suggests that, as input to Eq. (6), the surface pressure jump in Eq. (13) can 
be simplified to 

/ °0 

w{k 1 ,0)g(y 1 ,k 1 ,0)e' ,klUt dk 1 (18) 

-OO 

Note that the integration will be explicitly performed when Eq. (18) is input to Eq. (6). Only the 
A ?2 integration is neglected. Furthermore, the evaluation of the surface pressure terms at retarded 
time t — r/co will be executed in an exact fashion, as prescribed by Eq. (6). Note also that, as a 
function of observer time £, the surface pressure in Eq. (18) is constant in — b < < b. Although 

the present analysis involves an airfoil of finite span for which end effects are expected, Amiet [12] 
argues that these end effects are negligible except at very low frequencies. Therefore, the surface 
pressure in Eq. (18) will be applied to the entire airfoil surface, including the regions near the tips. 

The evaluation of the surface pressure in Eq. (18) is accomplished by first recognizing the 
turbulent fluctuations as a stochastic process. This process can be approximated by a truncated 
series whose coefficients are chosen such that the autocorrelation of the series forms a Fourier 
transform pair with its power spectrum (see, e.g., [24]). This property is achieved by evaluating 
the spectral amplitudes w(ki : fe) as a function of the power spectral density (PSD) of w{x\, X 2 , t). 
To this end, the infinite wave number domain, — oo < k\ < oo, in Eq. (18) is integrally discretized 
and truncated such that ki~N < ki^ n < fcqjv- The largest wave number ki^N represents an “upper 
cutoff” wave number, beyond which the spectral density amplitude is considered negligible or is 
out of range of experimental measurement. The unsteady surface pressure jump in Eq. (18) is then 
approximated by 

N 

& p (yi,t) « 2-kpqU ^2 A nfi e^ n g(yi,k ljn ,0)e lkl - nUt (19a) 

«=- N 

ki,n — Afci , n = 0, d=l, ±2, . . . , =L7V 
Aki = ki^/N 
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The phase angles {4^} are independent random variables uniformly distributed on [0, 27r]. The 
gust amplitudes {^, 0 } are related to the PSD of the turbulence and will be discussed momentarily. 
The transfer function g(yu 0) is the limiting case of Adamczyk’s [21] dual wave number transfer 
function and is given by 

g(yi,ki, n ,0) = ( <71 + <72 )e~ i [^(i-M) yi /c + ,/4-- kl ,n] (19b) 


where g\ and g 2, now functions of k\ jW , are otherwise identical to Eq. (9f), 


gi(yuh,n) 


1 

[2ir&i >n j/i(l + M)/C]% 


92 {yi,h,n) 


-l + (l+Q£*[4/Ul-yi/C)] 

[27rfc 1)n (l+M)]5 


(19c) 


and E* is the same complex combination of Fresnel integrals as in Eq. (9g). The low frequency 
transfer function is not used in these broadband predictions because the experimental facility is 
anechoic for frequencies above 200 Hz, and the parameter fi is greater than 0.4 at this frequency 
or above for all five tunnel speeds. 

The gust amplitudes {A n $ } are evaluated by 


— [ S ww (k\^ n , 0 ) Ak\ ] 


(19d) 


where S ww (ki,k 2 ) is interpreted as a three-dimensional PSD of w , into which all normal wave 
numbers are integrated, 

/ DO 

Sww (ki,k 2 ,h)dk 3 

-OO 

Although Ak 2 does not appear in Eq. (19d), the expression in Eq. (19a) still represents a dimen- 
sionally consistent surface pressure with respect to its role as input to the acoustic formulation in 
Eq. (6). This fact arises from the previous acoustic analysis in which k 2 was explicitly integrated in 
order to obtain the expression for the input surface pressure jump in Eq. (18). Moreover, the acous- 
tic analysis in Eqs. (14)-(17) alleviates the need for the scale factor that was previously reported 
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in [25]. It is important to note that the broadband predictions to be discussed in the following 
subsection are generated without any reliance on a scale factor or other arbitrary constant. 

As this grid generated turbulence is assumed to be homogeneous and isotropic, the PSD can be 
related to von Karman’s formula [20]. The formula for S ww (k\, fe) used in the present calculations 
is derived from von Karman’s energy spectrum in [13]. 


^ww (^1 •) ) — 


4 v? k\ + 

9^ fcf [l + jfe2 + fc2]| 


where u 2 is the streamwise mean-square turbulence, and 


l-h 12 k - 
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and C\ is the streamwise integral length scale 


Ci = 


POO 

I 

J 0 


dx i 


(19e) 


(19f) 


where R ww is the upwash correlation function defined by 


Rww{r) = {w(xi,t ) , w(xi +r,t) ) = lim ^ [ w(xi,t) w(xi +r,t) dt (19g) 

T-¥ oo T J o 

Note that the definition in Eq. (19g) implies that the stochastic process is assumed ergodic. 

The two-component formula in Eq. (19e) was derived in [13] by integrating the von Karman 
energy spectrum over all components. The &2 component is then set to zero for the present 
calculations. Values for u 2 and C\ are determined by experimental measurement. Fink [26] reports 
that the turbulence intensity that results from the grid used in this experiment can be approximated 
by the empirical formula 


U* ) 2 

U 


= 0.04 


U 

Urf 


I - 0.2 


(19h) 


where the reference speed is ?7 re f = 60 m/s. Fink also reports a measured value for the integral 
length scale to be C\ = 3.2 cm. 

All of the above expressions and measurements are incorporated into Eq. (19a). The final 
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representation for the unsteady broadband pressure on the airfoil surface is then given by the real 
part of —A P in Eq. (19a). This broadband surface pressure is used as input to Formulation IB to 
predict the acoustic pressure p f (x, t) at the microphone location. 

Using symmetry arguments and algebraic manipulation, the indicial bounds for the surface pres- 
sure’s spectral representation are altered so that the domain includes only positive wave numbers. 
As input to Eq. (6), the resulting real- valued surface pressure can be written 

N 

p(yi,r) = -AwpoU ^2 A nfi [B n cos(ki yn UT + - a n ) + D n sm(ki yn UT + - a n ) ] (20a) 

n—1 

where the upwash amplitudes A n $ are evaluated by Eq. (19d)-(19h), and the quantities B n , D n , 
and a n are given by 
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(20b) 


where gi(yu ki )U ) is evaluated in Eq. (19c), and C(^ n ) and <S(£ n ) are the Fresnel cosine and sine 
integrals in Eq. (9g) with £ n = 4// n (l — yi/C). The summation in Eq. (20a) begins at n = 1 
because 5^(0, 0) = 0, by Eq. (19e). 


4.3 Time Domain Predictions 

The lower frequency bound, and therefore the fundamental frequency, for all five calculations 
is chosen to be 10 Hz. The upper frequency for the predictions is chosen according to the upper 
frequency for which measurements are available for each tunnel speed. For U = 40, 60, and 90 m/s, 
the upper bound is /jv = 2.5 kHz. The upper bounds for U = 120 and 165 m/s are /at = 3.5 kHz 
and 4.5 kHz, respectively. For all five calculations, the numerical bandwidth is A / = 10 Hz. Each 
calculation is performed for one period of the lowest frequency, F = 0.1 s. The numerical solution 
is sampled at the Nyquist frequency, i.e., At = T/27V, where N is the upper cut-off index in Eq. 
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(20a). The calculation is performed on a 100 x 230 surface grid with the streamwise integration 
interval restricted as in Eq. (11), and with sufficient grid clustering near the leading and trailing 
edges to allow for e = 0.003. As in the constant frequency case, because the surface pressure is 
cast in only the one spatial variable 2 / 1 , and the observer location is symmetric relative to the 
airfoil span, the acoustic predictions were found to be relatively insensitive to the discretization in 
y 2 ‘ Again, the current discretization is employed as an example for a more general case, such as 
when the spanwise differences in retarded time become important for an observer whose location 
is asymmetric relative to the airfoil geometry. 

The experimental microphone position for which comparisons are made is at a distance of 2.25 
m from the model, and at an angle of 90 degrees relative to the geometric center of the model. 
The coordinate system for the calculation is such that the ^ 2 -axis coincides with the center span 
line, so that the microphone position is in the plane X 2 = 0, as in Fig. 7. The measured observer 
position for the prediction is, then, x — [0.115, 0, 2.25 ] T in meters. 

The position of the microphone relative to the airfoil is corrected for refraction caused by the 
presence of a shear layer that forms downstream of the upper lip of the square nozzle exit and is 
positioned between the model and the microphone (Fig. 9). This correction is based on geometrical 
acoustics with an assumption of a zero thickness shear layer, and is reported by Amiet [27]. Shear 
layer corrections that are based on such simple formulations are reasonable for the present case 
with the microphone directly above the source. The required correction in the microphone position 
is significant. At a measured angle of 90 degrees, the corrected angles ranged from approximately 
84.5 degrees for U = 40 m/s to 68.2 degrees for U = 165 m/s. 

In addition, the amplitude of the radiated noise is also corrected for the presence of the shear 
layer; however at a measured angle of 90 degrees the amplitude correction is not significant, espe- 
cially for the lower tunnel speeds. The computed sound pressures p f (x, t) were corrected by factors 
ranging from approximately 0.997 for U = 40 m/s to 0.942 for U = 165 m/s. The microphone 
position is corrected for the shear layer in a preprocessing step. After the far field noise is calculated 
at the corrected position, the results are then post-processed for amplitude correction. In this way, 
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the corrected predictions can be compared with the experimental results “as measured.” 

Fig. 10 shows the far field signal p f (x, t ) that is predicted by Formulation IB, for one fundamental 
period in time, at the experimental microphone location, for a tunnel speed of 165 m/s. Shear 
layer corrections for amplitude and directivity are included in this plot. In order to compare with 
experimental measurements in [8], the time domain results for all five tunnel speeds are Fourier 
analyzed and the resulting spectral amplitudes {Pn}n=i-> are used to compute the sound pressure 
level (SPL) spectrum of the far field radiation. These SPLs are determined by 


SPL(f n ) = 20 log 


Pref _ 


n = 1,2,. ..,7V 


( 21 ) 


where the reference pressure is P re f = 20 //Pa. The SPLs are converted to a 1 Hz bandwidth by 
reducing the values in Eq. (21) by 10 log(A /). This narrowband conversion is consistent with the 
experimental SPLs which were measured at a bandwidth of 55.7 Hz and reduced by 17.5 dB [8]. 

The predicted and measured far field SPLs for the five tunnel speeds are shown in Fig. 11. The 
experimental measurements in Fig. 11 were obtained by digitizing the measurements plotted in 
Figure 4 of Reference 8. The solid symbols represent those measurements for which the difference 
between the noise with and without the airfoil model was considered too small, and these mea- 
surements are therefore subject to greater uncertainty. The agreement with the measured data is 
very good, particularly in the range of 200 Hz to 1500 Hz, where the spectra are peaked. The rela- 
tively significant error in the prediction at the lowest tunnel speed is consistent with the low speed 
prediction in [8]. This error results from the fact that, as the air speed decreases, the gust wave 
lengths for the lower frequencies that dominate the noise will decrease and eventually become of 
the order of the airfoil thickness, thereby making the flat plate approximation increasingly invalid. 

The existence of multiple local extrema in a spectrum is often interpreted as the existence of 
multiple noise producing mechanisms; however, in the present case, the multiple extrema result 
from the changing placement of lobes, as a function of frequency, relative to a fixed observer. This 
phenomenon was described in Section 3.4 (Fig. 8). A similar trend in local extrema is also apparent 
in the experimental measurements at the two highest tunnel speeds, where the signal-to-noise ratio 
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was much higher. 

In Fig. 12, a comparison is made between the current time domain prediction and the frequency 
domain prediction of Paterson and Amiet [8], for the tunnel speed U = 165 m/s. The highest speed 
is chosen for the large frequency domain for which the spectrum is predicted. The frequency domain 
prediction in Fig. 12 was obtained by digitizing 51 points on the appropriate locus in Figure 4 of 
Reference 8. As the length of the square symbols is equivalent to 0.5 dB on the ordinate axis, clearly, 
the two predictions are very close throughout the entire frequency range, and nearly identical in 
the peak range of 200 Hz to 1500 Hz. That the current time domain predictions are so similar 
to the predictions in [8] is not surprising, as the current predictions rely upon the same unsteady 
surface pressure formulation as those in [8]. However, the decoupling of the aerodynamics and the 
acoustics that is apparent in Eq. (6) makes Formulation IB amenable to any type of input available, 
whether analytical, experimental, or computational. 

5. Statistical Formulation 

Often, when aeroacoustic experiments are performed, surface pressure correlations are extremely 
useful in characterizing noise due to an airflow over a model. Under certain conditions, Formulation 
IB is readily transformed into an expression that aids in the statistical analysis of broadband noise. 
Such an expression for the autocorrelation of the far field acoustic pressure is developed in this 
section. 

A far field autocorrelation formulation is derived for the case of a thin airfoil moving rectilinearly 
at constant velocity. The assumption of uniform velocity simplifies the algebra considerably, and 
the extension to general motion is straightforward. The following derivation follows the statistical 
analysis for jet noise presented by Morris and Farassat [28]. In the case of an observer in the 
acoustic far field, r A, the near field term (the second integral in Eq. (6)) can be ignored because 
of its 1 /r 2 proportionality. Also, for notational convenience, the spatial and temporal derivatives of 
pressure in the first integral are written as a single time derivative of pressure, p, when evaluated in 
a reference frame that is fixed relative to the medium at rest, as in Eq. (5b). In this fixed reference 
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frame, the far field acoustic pressure can be written 


47 xp 
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Eq. (22) is now transformed to the reference frame fixed to the airfoil, i.e. the wind tunnel reference 
frame. In this frame of reference, the emmision distance r ret has the following form: 


where 
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Therefore, in the wind tunnel reference frame, a time independent function K.(x,y) is defined by 
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Eq. (24) is an important result because it shows that, in the laboratory reference frame, the 
only time dependent quantities in Eq. (22) are p and p. Thus, Eq. (22) can be written in a form 
that shows the dependence of p f (x, t) on x , y, and t explicitly: 

4:irp'(x,t) = — [ JC(x,y)p(y,t-r ret /co)dy s - f M v {y) K.{x, y)p(y,t - r ret /co) dt (25) 
co Jf > o Jf= o 

where dys denotes an integral surface element dy \ dy-) . Using the notation (• , •) to represent an 
ensemble average, as in Eq. (19g), assuming ergodicity, the autocorrelation of the acoustic pressure 
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is evaluated as follows: 
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All primed quantities on the right-hand side of Eq. (26) indicate replacement of the variable y 
with y', e.g., K.' = K.(x,y'). 

Now, let y and y' be related by y' y -\ ff, where ff = [ 771 , 772 , 0] r is a vector in the plane of 
the airfoil. Furthermore, it is assumed that \ff\ <C r ret , in which case /) el can be approximated by 
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(27b) 


Two autocorrelation functions, Fm, and F pp , and a cross-correlation function F pp are therefore 


defined by 


Fpp{y 

o ff,f) 

= (p{$,t) 

,p(y + ff,T + f)) 


Fppii 

6 ff, f) 

= (p{y,T) 

,p(y + ff,T + f ) ) 

(28) 

Fpp{i 

6 ff,f) 

= (p{y,T) 

,i>(y + ff,T + f)) 



Using Eqs. (27) and (28), the ensemble-averaged quantities on the right-hand side of Eq. (26) have 
the following meaning: 
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(p(y, t - ^ret/co) ,p(y', t + f- r ret /co -ff-Vy r ret/co) ) = ^(y; 7?,f - 7? • r ret /co ) 

(p(y 5 1 — ^ret/co) ,p{y',t + t — r' ret /co — ff ■ r ret /c 0 ) ) = F pp (y; ff,f - ff • V y r r&t /c 0 ) (29) 

{p(y,t- net/co) + f - rJ et /co - ff- V y r ret /co) ) = F PP (y ; r?,f -rf- V y r ret /c 0 ) 

(p(i?,t ~ rret/co) ,p(y',t + f ~r' ret /co -ff ■ V y r ret /co)) = F pp (y ; -ff, -t + ff- V y r ret /co) 


Eq. (26) can thus be rewritten using the above correlation functions explicitly: 


16 7r 2 ( p'(x,t ) ,p'(x,t + r) ) = 

"4 [. dy s K.(x,y) [ 1C(x,y + ff) Ff,p(y-, ff,f - rf- V y r ret /c 0 ) dff 

c O J /> 0 «/ -A c 

- — [ dy s K(x, y) [ M v (y + ff) K{x, y + ff) Fpp(y ; -ff, -f + ff- V y r re t/co) di 

co y/>o Jlc 

~ — f dt M u (y)K{x,y) f K(x,y + ff) F pp (y-, ff,f - ff ■ V y r ret /co) dff 
Co Jf = O JAc 


(30) 


+ / d£ M v (y)K{x,y) / M v (y + ff) JC(x,y + ff) F pp (y; ff,r - ff • V y r ret /c 0 ) di 

lf = 0 Jf = 0 

Here, A c and L c are correlation area and correlation length, respectively. Clearly, performing the 
integration in the four inner integrals on the right-hand side of Eq. (30) will result in four functions 
of x. y. and f. Calling these function F\ to F 4 in order of their appearance, Eq. (30) can be 
rewritten 


(p'(x,t),p'(x,t + T)) = \ [ JC(x,y)[Fi(x,y,T)-coF 2 (x,y,T)] 

c o Jf> 0 


dys 


CO Jf= 


/ M v (y) JC (x, y) [ Fa (x, y,r) - cqF a (x, y, r) ] > 

J f = 0 


(31) 


The usefulness of the analysis presented above depends on the knowledge of two autocorrelation 
functions and one cross-correlation function, all of which are expected to be obtained from experi- 
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mental measurement. If the surface pressure fluctuations are assumed to be stationary in time, then 
only the correlation function F pp (y : ; 77 , f) needs to be measured, as the other two correlation func- 
tions can be derived from it [29]. It is expected that, of the three correlation functions, F pp (y\ rf, f) 
will dominate the acoustic pressure, and therefore, F\(y\ ff : f) will be the only significant contrib- 
utor to the right-hand side of Eq. (31). Eq. (30) indicates that, at a location y on the airfoil, only 
the turbulent pressure fluctuations over a correlation area near that point will contribute to the 
sound at the observer location x . Furthermore, the time dependence of the correlation functions is 
displaced by ff • V^r^t/co because of the retarded time effect. This displacement is dependent upon 
airfoil Mach number and the radiation direction x — y (the visual direction). Perhaps as important 
as its direct application to the prediction of far field noise from surface pressure measurements, Eq. 
(31) can be used for qualitative analysis of broadband and trailing edge noise, which is the subject 
of ongoing research. 
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Concluding Remarks 

The prediction of broadband noise from rotating machinery and airframes is currently the 
subject of intense research in aeroacoustics. The physics of broadband noise generation are well 
understood as the result of the pioneering research of Howe [14, 30, 31], Amiet and coworkers [3, 
4, 8, 13, 19, 32], and Brooks and coworkers [33, 34]. The previous work of these aeroacousticians, 
and many others, has clearly demonstrated that any successful broadband loading noise prediction 
requires an understanding of two physical processes: the character of the time-dependent surface 
pressure that provides the acoustic source, and the manner in which that source gives rise to an 
acoustic signal. 

The most prevalent types of unsteady surface pressure associated with broadband loading noise 
are those that arise from the interaction of the airfoil with a gust or the passage of eddies near 
the trailing edge. Obtaining this fluctuating pressure distribution analytically, numerically, or 
experimentally is itself a difficult problem. For this reason, past researchers have most often resorted 
to modeling the surface pressure, using guidance from experiments to aid in the development of 
these models. Today, high resolution surface pressure fluctuations can be obtained from turbulence 
simulations in realistic situations where the airfoil geometry and kinematics are accurately modeled. 
Therefore, the improvement of the acoustic radiation model becomes an important research topic. 
In the past, acoustic radiation models were most often developed for airfoils in uniform rectilinear 
motion. In addition, other restrictive assumptions, such as far field positioning of the observer, 
were often used to simplify the acoustic analysis. 

The present work develops a simple and general acoustic result in the time domain, based on the 
solution of the loading noise term of the Ffowcs Williams-Hawkings equation. This new solution, 
called Formulation IB, is, to date, the simplest analytical result for the prediction of loading noise 
and is suitable for statistical analysis of broadband noise for a surface in general motion. The 
new formulation has been validated with time domain calculations that predict the noise due to 
incident turbulence on a NACA 0012 airfoil. The time domain predictions are found to be in 
excellent agreement with the frequency domain predictions of Paterson and Amiet [8] as well as 
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with their experimental measurements over a significant range of tunnel speed. These results are, 
to the authors’ knowledge, the first successful broadband predictions in the time domain. 

Because of its relative simplicity, Formulation IB lends itself to a straightforward extension that 
provides an autocorrelation of the acoustic pressure in terms of two autocorrelation functions and 
one cross-correlation function of the surface pressure. Both numerical and experimental surface 
pressure distributions can be used in this result to yield information on broadband noise statistics. 

The authors advocate the use of time domain methods in the prediction of broadband noise. 
Because of the decoupling of the aerodynamics from the acoustics, the chief advantage of time 
domain methods is their potential for direct use of time-dependent surface pressure statistics from 
experiments or computer simulations in the prediction of broadband noise. 
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Figure 2 Differential surface element in Eq. (4c). 
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Figure 3 Schematic for the const ant-frequency loading noise problem in Section 3; an infinitely 
thin rectangluar surface in rectilinear motion. 
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Figure 4 Second level grid (20 x 80) for convergence study in Section 3.2. 
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Figure 6 Velocity scaling properties of Formulation IB. 
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Figure 9 Schematic for incident turbulence experiment in [8]. 
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Figure 10 Predicted far field signal, U = 165 m/s; microphone at 90°, 2.25 m above airfoil center. 
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Figure 11 Predicted and measured far field noise spectra; microphone at 90°, 2.25 m above airfoil 
center; experimental data reproduced from [8]; solid symbols denote low signal-to- noise ratio. 
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Figure 12 Comparison of predicted far field noise spectra, U = 165 m/s. 
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